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We investigate the nonlinear optical susceptibilities of Mott insulators with the dynam- 
ical mean field approximation. The two-photon absorption (TPA) and the third-harmonic 
generation (THG) spectra are calculated, and the classification by the types of coupling to 
external fields shows different behavior from conventional semiconductors. The direct transi- 
tion terms are predominant both in the TPA and THG spectra, and the importance of taking 
all types of interaction with the external field into account is illustrated in connection with 
the THG spectrum and dc Kerr effect. The dependences of the TPA and THG spectra on the 
Coulomb interaction indicate a scaling relation. We apply this relation to the quantitative 
evaluation and obtain results comparable to those of experiments. 
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1. Introduction 

Several nonlinear optical responses have been observed in Mott insulators; the two-photon 
absorption (TPA), 1,2 the third-harmonic generation (THG) 3,4 and the electroreflectance spec- 
troscopy. 5 A notable point is that quasi one-dimensional (ID) Mott insulators show large non- 
linear responses in these measurements, compared with those of conventional semiconductors. 
On the other hand magnitudes of nonlinear responses in two-dimensional (2D) systems are 
comparable to those of conventional semiconductors, and then the dimensionality dependence 
of nonlinear susceptibilities has also attracted attention in Mott insulators. 2,4 However this 
does not mean that the 2D system does not need an explanation, because the origins of the 
optical gap in the band insulators and Mott insulators are different from each other and the 
theory of conventional semiconductors does not apply to Mott insulators. There exists detailed 
comparison between experiments and theory in conventional semiconductors. 6 By contrast, 
optical nonlinearities in Mott insulators have not yet been understood to that level. 

In the previous paper we derived the general formulation of the nonlinear optical suscep- 
tibility based on Green's function, and applied this to a calculation of the TPA spectrum of 
antiferromagnetic insulators with the Hartree-Fock approximation. The dimensionality de- 
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pendence of nonlinear susceptibilities was investigated and a semiquantitative estimation was 
made there. This calculation fails to include the damping effect, and the divergence arises 
at the band edge. This makes a quantitative estimation difficult and it is done with the av- 
eraged spectrum. (The damping effect is also important due to the experimental fact that 
the response time of Mott insulators is very fast. 1 ) Other approaches on nonlinear optical 
responses have been made with use of the numerical diagonalization method on small-sized 
systems. 8 This calculation consists of the discrete levels and dipole moments between them, 
and requires an artificial damping term. Although the qualitative reproduction of the di- 
mensionality dependence is made with this method, even rough estimation of magnitudes of 
nonlinear susceptibilities is not attempted. 

In this paper we study nonlinear susceptibilities of Mott insulators with the dynamical 
mean field approximation on the basis of the general formulation developed in ref. 7 The 
damping effect is naturally included within this method. We calculate the THG spectrum 
and dc Kerr effect as well as the TPA spectrum. It is shown that the direct transition term 
predominates in the TPA and THG spectra. This is not the case in the dc Kerr effect, but 
all types of processes are important in the same degree to form the oscillating structure. The 
scaling relations of the optical responses are derived, and the linear and nonlinear responses are 
proportional to the inverse of the square and the fourth power of the energy gap, respectively. 
According to this relation we obtain quantitative results of the TPA and THG spectra, which is 
comparable to experiments in the case that the value of the Coulomb interaction is somewhat 
larger than that of the bandwidth. 

We present our formulation for calculation in §2, and the results are shown in §3. Several 
vertex corrections to the nonlinear susceptibilities are considered in Appendix. We set h = 
c = 1 and the electric charge e is not written explicitly. These are restored in quantitative 
calculations. 

2. Formulation 

Firstly we show how the Mott insulating state is described in our calculation. We apply 
the dynamical mean field approximation (DMFA) to the single-band Hubbard model, 



(tij is the transfer integral and U indicates the on-site Coulomb interaction.) We do not use 
the notation 'theory' which is usually used in the dynamical mean field theory (DMFT), but 
adopt 'approximation' because we do not consider the limit of the dimensionality d — > oo. 
This implies the following. In DMFT the effective single-site action, 




(1) 



<ij>a 




/ dr dr'^4(r)g - 1 (r-r')c CT (r') + ^ / drn^n^r) (2) 
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is derived in the large dimension limit, d — > oo. 9 (/3 = 1/T and T is the temperature.) In 
our case we use this effective action in arbitrary lattice systems. This means we neglect the 
higher-order terms of the transfer integral other than the first term of S e g . This is the reason 
why we use DMFA instead of DMFT. In this case we do not need to scale the transfer integral 
by the factor of 1/s/d. 

Other processes in the calculation are the same as in DMFT. The self-energy is calculated 
with S c g as the functional of Qq, S[^o]- The Weiss function Qq is calculated by the following 
relation, 

G l (e n )=le n+f i -G^(e n ) (3) 

and 



with Green's function, 



k 



G k {e n ) = -. ^-y (5) 

i£n ~ik + V- S(e„) 

(Here e n = 7rT(2n — 1) and n is integer.) These functions are self-consistently determined, 
and the chemical potential /x, jiQ is fixed by the condition, n CT [G] = ra CT [^o] = 1/2 (this sets 
the system to be half- filled) . We make another approximation to solve S e ff- We calculate the 
self-energy within the second order perturbation, 10 

E(e„) = -U 2 T 2 Go{e n ,)g {e n , + uji)g (e n - wj). (6) 

n',l 

We use the following dispersion relation, 

= — 2i(cosfe r + i]cosk y ) + 4t'r]cosk x cosk y . (7) 

In numerical calculations below we put t = 1 and fix the next-nearest-neighbor hopping 
t' = 0.2 (results do not change if we vary t' moderately). We vary r\ as the dimensionality 
parameter from the 2D r] = 1.0 to the quasi ID r] = 0.1. 

Next we present the formulation of the nonlinear optical response functions. The third- 
order nonlinear susceptibility is determined by, 

X l3) {^,u 1 ,u 2 ) = 1 ; . (8) 

(The definitions of x (3) and K^> are given in ref. 7 ) Here, u> = uj\ + u>2 + <^3 , ^' = <^2 + <^3 
and oj" = UJ3. ui, uj 2 and ujs are frequencies of the external fields and take different values 
depending on various methods of measurements. K^> is classified by the types of the coupling 
to the external fields as follows, 

( U ,u/, u") = K%> + K%> + K%> . (9) 
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Each term is written as, 



K %> = lX I ^4Y.^GfGfT b + G^GfT 3 Gt + G^T l GfGt + T a GtGfGt), (10) 

k J <i,j> 

K %> =| E / wm^51^ G ? T » + G ^ + T ^ G ^ + E^Gfn + g«t 3 g a + T aG f G A ) 

k i j 

+ £ (GfGfT b + G?T jG A + TiGfGf)], 



<i,3> 



and 



h i j 



+ ^v k (GfT b + T a G A - GfT h - T b G A )\. 



(11) 



(12) 



Here G^ A = G*' A (e x ) (R and A mean the retarded and advanced, respectively), T x = 
t&nh(e x /2T)ImG^(e x ) and = d^/dk. e x = e + lo x and uj a = u\ + uj2 + <^3, = 0, 
uji = io\ + ll>2, uj\ + CJ3 or 0^2+^3, = wi, LU2 or uj%. The diagrammatic representations are 
given in Fig. 1 of ref. 7 ; Fig. 1(a), (b) and (c,d,e) for K^ 4> , K^ 3> and K^ 2> , respectively. 
In this formulation vertex corrections are omitted, and these are discussed in Appendix. 

3. Results 

3. 1 The analysis of spectrum 

The numerical results shown below are calculated with eqs. (10,11,12). The vertex correc- 
tions are not included, which are small compared to these terms as indicated in Appendix. 

The decomposition of ImK^ to ImK^j A> , ImK^ 3> and Imif^ 2> in the case of the 
TPA spectrum {uj\ = —LO2 = uj^ = uj) is shown in Fig. 1. (We fix the temperature T = 0.036 
hereafter, and this parameter is not considered to be important because of u,U,t S> T.) The 

(3) (3) (3) 

predominance of K y < - 2> over ^<j4> an d -^<j3> ^ s peculiar to Mott insulators, in contrast 

with conventional semiconductors where K^ 2> vanishes except for the self-transition. 11 ' 12 

The existence of K^ 2> i n the- TPA spectrum depends on the origin of the gap, 7 and the 

difference in magnitude of three ImK^ is understood by writing expressions explicitly as 

follows. 

Imi4 3 j 4> (u;,0,u;) = \ J ^^( tanh ^^ ~ tanh l^) ^Reh-^Re 

k 

+ ^tanh-^ tanh— ^ I e+ul I e {I e+ul I e — R t+UJ R t — R e+ 2 W R e +u — ReRe-u)]- 

(13) 
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Fig. 1. The decomposition of Imi^ 3 ) in the case of the TPA spectrum. U = 12 and -q = 0.2 'sum' 
indicates the sum of three terms. 



In^3>(">0.") = |£ / 



k v\ [ ( tanh 



3^ J 2vr dk ftL V 



e — uj e + uj 

tanh — ) Ie+uRele-u 



2T 



2T 



+ ^tanh^ - tanh^^ I e+LU I e (R e+UJ + R e + R e +2uj/2 + R e 



(14) 



V2)]- 



^>(-0,.)4E/|[(^) 2 (tann 



e-uj e + uj 

-TfjT ~ tanner" ) I €+W I € - L 



d 2 Vk ( . e 
+2 W fc l tanh 2T 



(15) 



Here J e = ImGf(e) and J? £ = ReGf (e). We consider the case of uj ~ Z7/2 in the TPA 
spectrum. These expressions indicate that the second terms of these three equations are small 
due to the factor I e +uh = I m G^(e + w)ImG^(e). (If one of lmG R takes large values, the 
other has small values owing to the absence of the spectrum.) Then we consider the first 
terms in these expressions. The existence of R e is the reason for the smallness of K^ 4> 

(3) (3) 

and K^j 3> compared with the direct transition term K y < j 2> . The former two cases includes 
virtually excited states in the optical process, and R e expresses this excitation. Although 
= ImGf(e + w)ImG|?(e — uj) can take large values around e ~ 0, R e is roughly 
proportional to 1/U in this region and is small. This explains results of Fig. 1. 

The decomposition of to K^ 4> , K^ 3> and K^ 2> in the case of the THG spectrum 

(3) (3) (3) 

(uji = L02 = oj 3 = uj) is shown in Fig. 2. The predominance of -fQj 2 > over 

K% 4> and K^j 

3> 1S 

the same as the case of the TPA spectrum, and the reason for this is also the same. (Here we 
consider the case of uj ~ U/3.) If we write the expressions of explicitly, we can find that 
the factor like ImGf (e+3w)ImGf (e) exists in K^ 2> . Then K^ 2> takes larger values than the 
other two terms, which include the nonresonant R e term. In the THG spectrum the existence 
of the real part Rex*- 3 -* makes it inevitable to calculate all three terms of consistently, 
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Fig. 2. The decomposition of (a) the real and (b) imaginary part of in the case of the THG 
spectrum. U = 12 and rj = 0.2. 'sum' indicates the sum of three terms. 



0.0014 



>5 




Fig. 3. The decomposition of IXthgI m tne 

case of the THG spectrum. |x<j4>L Ix <j3> \ and |x<,2>l 



are calculated with K^ 4> , K^j 3> and K K ^ 2> , respectively, 'all' means |x<j4> + X<j3> + X<j2>l- 
U = 12 and rj = 0.2. 



,(3) 



,(3) 



(3) 



.(3) 



X3) 



especially for small uj. If we calculate IxthgI om y with K<j4>i ^<%> or ^<j2> separately, 
each of Ix<j4 j3 j2>l diverges at small uj as shown in Fig. 3. The cancellation among three 



K 



(3) 

<j4,j3j'2> 



occurs at small uj, and we obtain convergence only if the summation of these three 



terms is taken. (This cancellation is the nonlinear analogue of that between the paramagnetic 
and diamagnetic terms in the linear response. It is unaffected by vertex corrections owing to 
the absence of the momentum-dependence in the self-energy.) This shows the importance of 
taking all three terms into account. The convergent behavior is related to that of the Drude 
weight, which is defined as D := 7ruTmcr| w ^o (a is the conductivity) and D = for T — > 
in insulators. 13 The nonlinear correction is written as Iniff' 3 ' = —HeK^/uj 3 . Therefore the 
nonlinear correction to the Drude weight would be divergent if KeK^ took finite values. 
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Fig. 4. The decomposition of ImK *- 3 '(w, 0, — Auj). Aui = 0.05, U = 12 and n = 0.2. 'sum' indicates 
the sum of three terms and the inset shows this result separately because of the difference in scales. 



(Strictly speaking, the term which is proportional to ex.p(—E g /T) (E g is the energy gap) 
remains as in the linear response, but this is vanishingly small for E g 3> T.) 

The decomposition of ImK^ to ImK^j 4> , Imi^^- 3> and Imi^^- 2> in an approximate 
case of the dc Kerr effect (uj\ = uj, UJ2 = —<^3 = Acj) is shown in Fig. 4. (It should be 
Auj — ► 0, but we apply the finite difference to = —K^/(uj 2 Auj 2 ).) In contrast to the 
above two cases, all of i4 3 j4» K <j3> and K <j2> contributes to in the same degree. The 
reason for this is that ReG^(e) does not necessarily locate at the nonresonant state, which 
is understood by writing the set of frequencies; (uJi,LUj) = (0, Au), (0, — Auj), (uj + Auj, Auj), 
{uj + Auj,uj), (uj — Auj,— Auj), (uj — Auj,uj). As shown in the inset the summation of these 
three terms is smaller than each of them by two orders of magnitude. All these terms are 
required to reproduce the characteristic oscillating structure similar to that observed in the 
electroreflectance spectroscopy. 

It is known that sum rules hold in the nonlinear response. 14 ' 15 The relation, 
J °° uje2 L (uJ, —uj', uj')duj = holds for the TPA spectrum. Here, e^ L (w, —uj', uj') is the imaginary 
part of the complex dielectric function. If we treat the above three terms of separately, we 
will violate this relation. The appearance of the oscillating structure in the dc Kerr effect as 
shown above is another example of the necessity to consider all these terms in (uj' = Auj 
in this case). A previous calculation do not take these terms into account properly. 16 They 

(3) 

neglect the predominant term K y K - 2> , and also fail to treat the divergence at small frequency 

(3) 

region carefully in a calculation of the real part of Xthg 

3.2 The dependences of nonlinear susceptibilities on the Coulomb interaction 

We show the dependences of the nonlinear susceptibilities on U/t, U/W and r/. (Here 
W is the bare bandwidth and is a function of t' and n.) The dependence of the integral 
of the linear absorption spectrum (a = J ujlmx^ (uj)duj) on t/U, W/U with several values 
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t/U 

Fig. 5. The dependence of a = J ujlmx^ (uj)duj on t/U. The inset shows the dependence of the same 
quantities on W/U. 



of r\ is shown in Fig. 5. (The value of U at which the Mott transition occurs depends on 
77, and these are U ~ 9.1,9.2,10.2,13.1 for 77 = 0.1,0.2,0.5,1.0, respectively.) The relation 
a oc 1/U holds, which is consistent with the sum rule for the linear absorption. 1 If we put 
the lattice constant a = 5 Aand U = 2 eV (The reason why we take this value is that the 
linear absorption spectrum in experiments peaks around this energy and in our simple model 
the spectrum always has the peak around U.), we get a^Leak — 0.98,0.76 x 10 5 cm -1 at 



U = 13.5 for 77 = 0.1, 1.0, respectively (here a 



Mi 



UJ 



A-Kojlmx^ (w) / c and c is the velocity of 



light which is written explicitly for the quantitative estimation) . These are almost comparable 
to the results of experiments which are | pca k — 4, 1 x 10 5 cm^ 1 in quasi ID and 2D systems, 
respectively. 2 The relation ulmx^ oc 1/U indicates that we expect a moderate enhancement 
of c/ 1 ) for smaller U. 

(3) 

The dependences of the peak of the TPA spectrum (lm.x T p A multiplied by to) on t/U and 
W/U are shown in Fig. 6. The relation wImx^p A oc 1/U 3 holds approximately. (It deviates 
slightly from 1/U 3 for smaller U, and the results are rather proportional to 1//7" 3 - 5 . This 
is because the peaks of the TPA spectrum shift to lower energies.) If we put the lattice 
constant a = 5 Aand U = 2 eV, we get Imx§L A ^ 0.0155, 0.0133 x 10~ 9 esu at U = 13.5 for 
77 = 0.1, 1.0, respectively. If we extrapolate the relation Imx^p A oc 1/U A for smaller U, we will 
obtain Imx^pA — 1-0,0.1 x 10~ 9 esu at U = 4.76, 77 = 0.1 and U = 8.15, 77 = 1.0, respectively. 

The dependences of the peak of the THG spectrum on t/U and W/U are shown in Fig. 7. 
The relation IxthgI °^ l/U^ holds. If we set parameters same as above to evaluate IxthgI 
quantitatively, we get IxthpI - 0.0217,0.0162 x 10~ 9 esu at U = 13.5 for 77 = 0.1,1.0, 
re5 p e ctive ly . If we a™ that the r ela«on | X « G I « W hohds tor am a!le r U, we w,U ohtain 
IxthgI — 1-0,0.1 x 10~ 9 esu, at U = 5.18, 77 = 0.1 and U = 8.56, rj = 1.0, respectively. 

These results indicate that the dependence of the susceptibility on 77 is rather weak, at 
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Fig. 7. The dependences of |xTHGlp cak 011 t/U- The inset shows the dependence of the same quantities 
on W/U. 

least with t/U fixed. On the other hand it is strongly dependent on r\ in the case of W/U fixed, 
and this is because the bandwidth W is a function of rj. The experimental results indicate that 
X (3) n 1.0,0.1 x 10" 9 esu for quasi ID and 2D systems, respectively. 2 ' 4 Our calculation shows 
that it is possible to obtain comparable to those of experiments in the case of U > W 
(actually W = 4.4 and 8.0 for r\ = 0.1 and 1.0, respectively). However this is based on the 
condition that we can extrapolate scaling relations for smaller U, and we discuss this point 
in §4. We find that the dependences of x on t' is weak with the moderate variation of t' . 

In experiments the nonlinear susceptibility in the quasi ID system is one order of magni- 
tude larger than that in the 2D system. Our result does not show so much difference between 
7] = 0.1 and 7] = 1.0 with fixed t/U. Although the improvement should be done on DMFA 
especially in quasi ID systems, this is partly explained by the behavior of the density of states, 
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Fig. 8. The density of states p(e) = — J2k I m( ^fc ( e )/ 7r with several values of U and rj. 
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Fig. 9. (a) Imx^p A and (b) \Xthg\ wn ^ n scvera l values of U and rj. 



which is shown in Fig. 8. The experiment on the linear absorption spectrum indicates that the 
band-edges of the spectrum are almost same in both systems. This means that the nonlinear 
susceptibilities to be compared should have the same band-edge in the density of states. There- 
fore we compare the nonlinear susceptibilities at U = 10.5, rj = 0.1 and U = 13.5, rj = 1.0 as an 
example having such properties. The TPA and THG spectra are shown in Fig. 9. According 
the scaling relation Xtpa thg ^ V^ 4 ' a slight change of U brings about large variations in the 
nonlinear optical susceptibilities. On the other hand the linear absorption spectrum does not 
change considerably because of x*- 1 -* °c 1/U 2 . Consequently the ratio of |?j=o.i to x^|r?=i .0 
becomes much larger than that of x( 1 )|^ =0 .i to x^|tj=i.0j which resembles the observations in 
experiments. 

The scaling relation in semiconductors shows that Xtpa ^ l/^- 6,12 (Eg is the energy 
gap.) Although this is similar to our result, this does not mean that both Mott insulators 
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and conventional semiconductors obey the same scaling relation because the dominant terms 
in are different between these materials as mentioned in §3.1. In spite of this fact, the 
difference in the magnitude of the nonlinear susceptibility between these materials is partly 
explained as follows. For the low dimensional systems the gap edge of the density of states 
is steeper than that of more high dimensional systems as shown in Fig. 8. This enhances the 
magnitude of the optical susceptibility in quasi ID systems, compared to that conventional 
semiconductors. 

4. Summary and Discussion 

We calculate nonlinear optical susceptibilities with DMFA on the basis of the general for- 
mulation of nonlinear response developed in a previous paper. The direct transition term is 
predominant in the TPA and THG spectra, which is contrary to conventional semiconductors. 
This is because the transition to the nonresonant intermediate states gives small contribution 
to due to the strong correlation. On the other hand the origin of the band gap in semicon- 
ductors makes the direct transition negligible in x^ ■ In spite of these facts our result shows 
that as a function of the energy gap the scaling relation in Mott insulators behaves similarly 
as that of conventional semiconductors. A semiquantitative evaluation of nonlinear suscepti- 
bilities is carried out and shows that results are comparable to those of experiments on the 
condition that the value of the Coulomb interaction is somewhat larger than the bandwidth. 
The magnitude of Imx^pA an d IXthgI takes similar values with each other, which is also 
indicated by experiments. These are not clarified in previous works for small systems which 
are diagonalized numerically. The scaling relation based on DMFA also shows that the smaller 
U is favorable to the larger x® as in the Hartree-Fock calculation, which is contrary to the 
scenario of a large optical nonlinearity based on the spin-charge separation. 8 (The spin-charge 
separation holds approximately and is preferred at large U/t. The validity of the spin-charge 
separation as an explanation for the large optical nonlinearity can be judged partly from the 
dependence of x*- 3 -* on parameters like U/t.) 

One of our conclusions is dependent on the assumption that the scaling relation holds for 
smaller U. Here we discuss on this point and a possible modification. The main reason why 
the relations x^ oc 1/U 2 and Xtpa thg ^ V^ 4 hold is as follows. By definition x^ oc 1/oj 2 
and Xtpa thg ^ V^ 4 - This leads to the above [/-dependences on the condition that the 
[/-dependences of K^> and ^tpathg are wea ^ an( i the optical gap scales with U. The 
calculation here indicates that this property holds at least within our approximation. However 
there is some room for improvement with respect to the description of the Mott insulator. The 
Hubbard model is considered to have the Mott transition at smaller values of U than those 
of a calculation presented here. This is the case especially in the model with r\ = 0.0, which is 
the ID system and should be the Mott insulator even as U — > 0. 18 The improvement should 
be done on this point to examine the dimensionality dependences and the scaling relation for 
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smaller U (for example, an expansion to include /c-dependence of the self-energy 19 ). 
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Appendix: Vertex corrections 

The correction to vertices v k and d 2 Vk/dk 2 vanishes as in ref. 20 ' 21 owing to the inversion 
symmetry. On the other hand it is not known to what extent the correction to vertices dv k /dk 
and v\ contributes to x^ ■ The diagrams and equations of this type of vertices are similar to 
those of Fig. 1 (f,g,h) and §3.2 in ref. 7 The vertex correction to the predominant term in 
is written as, 

K^(m) = -T 2 £ E Gfc ( e « + "l)^G k (e n )T { e n , e' n , W ,) £ G k ,{e n> + U i)-^-G kl (e n >). 

n,n' k k' 

(A-l) 

Here r(e n , t n <\ loi) is the reducible four-point vertex. If we consider the second-order 
perturbation term as an irreducible four-point vertex /(6 n , e n i it is written as 

I(e n ,e n r,ui) = 2x(e n - e n >) + cf>(e n - e n >) (x{^i) = -U 2 TJ2 k)n G k (e n + ui)G k (e n ) and 
4>(uii) = —U 2 T'}2 kn G k {uJi — e n )G k {e n )). From the expression we anticipate that the ver- 
tex correction is small in the case that the dependence of on frequency is weak. It is 
because the particle-hole symmetry holds approximately. This can be verified by the numeri- 
cal calculation which shows that the vertex correction is smaller than ordinary terms by two 
orders of magnitude. 

In contrast to this, the nearest-neighbor interaction (V) is considered to be important in 
optical responses because the excitons can be formed by the final-states interaction. Therefore 
we consider the vertex correction by the nearest-neighbor interaction. The formulation is 
similar to that of §3.2 in ref. 7 and we consider only the Fock term with this interaction. The 

(3) (3) 

vertex correction to the predominant term ImK^' 2> , lmK y < j 2> itself and the summation of 
both terms of the TPA spectrum are shown in Fig. A-l. The vertex correction shifts the 
spectrum to lower energy. This effect is rather small compared to that of antiferromagnetic 
insulators with the Hartree-Fock approximation because the damping effect is included in 
DMFA. Although the value of V/U is not known (V/U ~ 0.22 in Fig. A-l is considered to 
be a large value), this type of the vertex correction will increase the values of x^ of §3.2 in 
some degree. 
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Fig. A-l. The vertex correction to Kj,p A by the nearest-neighbor Coulomb interaction V with the 
Fock approximation. U = 13.5 and V = 3.0. 'vc' and 'no vc' means the vertex correction term 
and ImK^, A without the vertex correction, and 'sum' means the summation of both terms. The 
vertical and horizontal axes are scaled with three times and twice values, respectively. 
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